Abstract 



We study the canonical structure of the SU(N) non-linear cr-model in a 
polynomial, first-order representation. The fundamental variables in this de- 
scription are a non-Abelian vector field and a non-Abelian antisymmetric 
tensor field 9^, which constrains L M to be a 'pure gauge' (F^ U (L) = 0) field. 
The second-class constraints that appear as a consequence of the first-order 
nature of the Lagrangian are solved, and the reduced phase-space variables 
explicitly found. We also treat the first-class constraints due to the gauge- 
invariance under transformations of the antisymmetric tensor field, construct- 
ing the corresponding most general gauge-invariant functionals, which are 
used to describe the dynamics of the physical degrees of freedom. We present 
these results in 1 + 1, 2 + 1 and 3 + 1 dimensions, mentioning some properties 
of the d + 1-dimensional case. We show that there is a kind of duality between 
this description of the non-linear u-model and the massless Yang-Mills theory. 
This duality is further extended to more general first-class systems. 
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1 Introduction 



One of the distinctive properties of the non-linear cr-model 0, is that its dynamical 
variables belong to a non-linear manifold thus realizing the symmetry group in 
a non-linear fashion ||. Whence, either the Lagrangian becomes non-polynomial in 
terms of unconstrained variables, or it becomes polynomial but in variables which 
satisfy a non-linear constraint. It is often convenient to work in a polynomial or 
'linearized' representation of the model. By this we mean an equivalent description 
where the symmetry is linearly realized, although the transformations now act on 
a different representation space. There is usually more than one way to construct 
such linearized representations. For example, in the O(N) models, where the field 
is an ^-component vector constrained to have constant modulus, a polynomial 
representation is constructed simply by introducing a Lagrange multiplier for that 
quadratic constraint. However, this simplicity is not present in general because the 
constraints required to define the manifold can be much more complex, like in the 
SU(N) groups. 

In references 0, |5j a polynomial representation of the non-linear cr-model was 
introduced; let us briefly explain it for the SU(N) model in d+1 dimensions. 

The usual presentation || of this model is in terms of an SU (N) field U(x), with 
Lagrangian density 

£ = ^g d -Hr{d^d^U) , (1) 

where g is a coupling constant with dimensions of mass (the constant f n in its 
application to Chiral Perturbation Theory in 3+1 dimensions). The polynomial 
description |5[] of this model is constructed in terms of a non-Abelian (SU (N)) vector 
field plus a non-Abelian antisymmetric tensor field 9^J^ with the Lagrangian 

1 To avoid the proliferation of indices, we frequently work in terms of the dual of 0^, which in 
1 + 1 is a pseudo-scalar, in 2 + 1 a pseudovector, etcetera. 
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density 

C = l -g 2 L,-L» + g6, v -F^{L) (2) 

where the fields L M and 6^ v are defined by their components in the basis of generators 
of the adjoint representation of the Lie algebra of SU(N); i.e., L^(x) is a vector with 
components L a ,a = 1, ■ ■ ■ , N 2 — 1, and analogously for 9^ u . The components of F^ u 
in the same basis are: F° V {L) = d^L a u -d u L°+g^ j ahc L\L c v . The dots mean SU(N) 
scalar product, for example: ■ = X^f 1 L a ^L^. The rf-dependent exponents in 
the factors of g are chosen in order to make the fields have the appropriate canonical 
dimension for each d. 

The Lagrange multiplier 9^ u imposes the constraint F^ U (L) = 0, which is equiv- 
alent || to = g~^UdfjlJ\ where U is an element of SU(N). When this is 
substituted back in (|2|), (P is obtained 0. This polynomial formulation could be 
thought of as a concrete Lagrangian realization of the Sugawara theory of cur- 
rents J7|, where all the dynamics is defined by the currents, the energy-momentum 
tensor, and their algebra. Indeed, corresponds to one of the conserved currents 
of the non-polynomial formulation, due to the invariance of £ under global (left) 
SU(N) transformations of U(x). The energy-momentum tensor for ([|) is indeed a 
function of only: 

T i» = g 2( L n . L v _ ~ g ^L 2 ) . (3) 

One can easily relate amplitudes with external legs of the field to the corre- 
sponding pions' scattering matrix elements, as shown in ref . || . It is also possible to 
relate off-shell Green's functions of the field U to the ones of the field L M , although 
this relation is non-local. As 

Lfx = g^Udjfl => D^U = 0, = 8^ + g'^L^ , (4) 

2 For a complete derivation of the equivalence between the theories defined by (|l|) and (Q) within 
the path integral framework, see ref. Q. 
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then U can be obtained at the point x by parallel transporting its value at spatial 
infinity, which we fix to be equal to the unit matrix^: 

U(x) = Vexp[-g^ [ dy^L^y)], (5) 

where V is the path-ordering operator ||, and the line-integral in the exponent is 
along a curve C x , a regular path starting at spatial infinity, and ending at x. The 
condition F^ u = guarantees that U is in fact invariant under deformations of C x 
which leave its endpoints unchanged. We can also construct products of two or more 
fields in a similar way, for example 

U(x 2 )U~ 1 (x 1 ) = Vexpi-g'^ f dy^(y)) , (6) 

where C x -> y is a continuous path from Xl to x 2 - This shows how [/-field correlation 
functions can in principle be calculated using Lagrangian (□); one has to evaluate, 
for example, the Wilson line (|J) in the theory defined by @. 

The classical equations of motion for the Lagrangian @ are 

L v {x) = -D^(x), 
F» V {L) = 0. (7) 

Taking the covariant divergence on both sides of the first equation of motion, and 
using the second one, one gets 

d-L(x) = . (8) 

Inserting L M = g^Ud^W in it yields the equations of motion for the usual 
non-polynomial Lagrangian ([IJ). Note that the solutions of (^) will, in general, 
contain arbitrary functions of the time. If we know a solution, performing on it the 
transformation: 

U^u(x) = £>'uv,(z) , (9) 

3 We identify (as usual) all the points at spatial infinity. 
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where u Pflu (x) is an arbitrary completely antisymmetric tensor field, will produce 
another solution, because D^ 1 6 UJ 9 fll/ vanishes as a consequence of the Bianchi identity 
for L^. Obviously d must be larger than one in order to this transformation be well 
defined, since at least three different indices are needed to have a Bianchi identity. 
This degeneracy in the equations of motion is due to the gauge invariance of the 
action under the transformations (Q). 

This gauge-invariance makes the quantization of the model interesting, and it 
will allow us to discuss some properties of the non-linear u-model from the (unusual) 
point of view of gauge systems. The Hamiltonian formulation of the model possesses 
a rich structure, since there are second-class constraints (£ is first-order), first-class 
constraints (for d > 1), and moreover they are reducible for d > 2. 

The structure of the paper is as follows: In section 2 we discuss the Hamiltonian 
formulation of the 1 + 1, 2 + 1 and 3 + 1 models, following the Dirac algorithm |J. In 
section 3 we construct the general gauge invariant functionals for the transformations 
generated by the first-class constraints found in section 2, and in section 4 we apply 
the Dirac's brackets method to the second-class system formed by the first-class 
constraints plus some canonical gauge-fixing conditions. In section 5 we present our 
conclusions. 

In Appendix we discuss a duality relationship between first-class systems, which 
generalizes a property we discuss for the 2 + 1-dimensional model. 

2 Hamiltonian formalism and constraints 

4 This kind of symmetry also appears when considering the dynamics of a two-form gauge field, 
see for example references O] 
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2.1 1 + 1 dimensions 



From Section 1, the polynomial Lagrangian in 1 + 1 dimensions becomes 

C = l - g ^L^+ l -g9e^{L) (10) 

where 9 is a pseudoscalar field. It is evident that there is no gauge symmetry in this 
case. Thus there will not be first-class constraints in the Hamiltonian formulation. 
However, there are second-class constraints, because L is of first-order in the deriva- 
tives. This property will also appear in higher dimensions, so we will only discuss it 
in some detail for this case. To start with, we rewrite (O) in a more explicit form 

C = l -g 2 L a L a , - \g>L\Ll + g6 a d L a x - gd a d x L% + tfPf**L\L\ . (11) 

Next we define the canonical momenta, where the primary constraints appear: 

* s mm K9r(x) 
™ s am" (12) 

and the canonical Hamiltonian becomes 

H = J dx[~g 2 L a L a + \g 2 L\L\ - gL^Of] , (13) 

where 

(D^y = d 1 e a + g f abc L\e c . (14) 



The 'total' Hamiltonian is constructed as usual, adding to ([131 ) a Lagrange multi- 



plier term for each of the primary constraints (|12l). Following the Dirac's algorithm 
one more constraint is obtained: 

gLl{x) « -{D x 0) a (x) (15) 



and the Lagrange multipliers become fully determined. The full set of (primary plus 
secondary) constraints is second-class, and its particular form allows us to eliminate 
the canonical pairs of Lq and 9 a , thus effectively eliminating the associated degrees of 
freedom. The Dirac bracket becomes equal to the Poisson bracket for the remaining 
degrees of freedom. The resulting Hamiltonian is 

H = J dx^Dnn • Ditti + ^g 2 L x ■ L x \ , (16) 

with canonical brackets between the L°'s and their momenta 7r". Thus these two 
variables become symplectic coordinates on the reduced phase-space, or constraint 
surface. 

2.2 2 + 1 dimensions 

The polynomial Lagrangian in this case becomes, 

£ = ±g 2 L„ .L»+ l -g6,- e^ x F uX (L). (17) 

The constraint algorithmfj produces the 2 + 1 analogous of the second-class con- 
straints we showed in Section 1, allowing us to eliminate the 0-component of L M and 
all the components of 6* M . However, there will remain a set of first-class constraints 

G a (x) = l -e jk F« k « , (18) 

with the first-class Hamiltonian 

H = J d 2 x[ 1 ^D J Ti j ■ D k 7i k + l -g 2 L 3 ■ L 3 \ . (19) 

They satisfy the algebra 

{G a (x),G b (y)} = 

{H,G a {x)} = V ab G\x) 
V ab = g-ir\D^ ) c {x). (20) 

5 The full details of the application of the Dirac algorithm to this system will be presented 
elsewhere. 
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Now we show in what sense we can relate the massless Yang-Mills theory to the 
non-linear cr-model in this formulation. The SU(N) Yang-Mills theory is defined by 
the Lagrangian 

Cym = -^(L)-F^(L), (21) 
which in the temporal gauge gives rise to the canonical Hamiltonian 

H = j d 2 x[^n 3 ■ nj + ~F jk ■ Fjk] , (22) 

and the first-class constraints ('Gauss' laws'): 

H a (x) = (D^) a (x) « , (23) 

which satisfy the SU(N) algebra 

{H a {x), H b (y)} = 8(x - y)f abc H c (x) . (24) 

Note that (^) can be rewritten as 

H = J d 2 x\^K 3 ■ -Kj + ^G-G], (25) 

where the G a 's are the ones defined in fli"8|) . We then see that the first-class systems 
corresponding to the Yang-Mills model and the non-linear cr-model can be related 
by: 1) Interchanging the constraints: 

H a (x) <-> G a (x) , (26) 

and 2) Interchanging Lj by ^iij in the non-derivative terms in the Hamiltonian. The 
generalization of this mapping is constructed in Appendix. 

2.3 3+1 dimensions 

After eliminating the second-class constraints, one obtains a first-class Hamiltonian 
which looks exactly like the one of the 2 + 1-dimensional case: 
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and the set of first-class constraints 

Gftx) = \e m F^{x) « . (28) 

Although the system seems to be the obvious generalization of the 2 + 1-dimensional 
one, there is an essential difference: The constraints (p8| ) are not all independent, 
but verify the Bianchi identity: 

(DjGjYix) = , Va. (29) 

This implies that the set of constraints is reducible, containing only two independent 
functions. The counting of degrees of freedom then gives 1 for the number of physical 
dynamical variables (1 = 3 — 2). We mention that the elimination of the second- 
class constraints applies in a similar way to the general d-dimensional case, and that 
the Hamiltonian and constraints are the obvious generalizations of ( p7|) and (|28D , 
respectively. Due to the existence of the Bianchi identity in general, the number of 
independent constraints in an arbitrary dimension is just enough to kill d — 1 out 
of the d degrees of freedom in H, leading to only one physical variable, as it should 
be for a model which describes the dynamics of a scalar field. 

3 Gauge invariant functionals 

Gauge invariant functionals^ are important from both the classical and quantum 
mechanical points of view. Classically, a complete set of gauge invariant functionals 
and their equations of motion completely determines the dynamics of the observable, 
i.e., physical degrees of freedom. In Quantum Mechanics, Dirac's method for first- 
class constraints defines the 'physical' subspace of the complete Hilbert space as the 
one whose state vectors are annihilated by the first-class constraints, i.e., the gauge- 

6 We assume the denomination 'gauge-invariant' to mean on-shell gauge invariance, i.e., the 
gauge-invariant functionals are invariant on the constraint surface. 
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invariant ones. In the Schroedinger representation, the physical subspace consists 
of gauge invariant functionals of the fields. 

To construct the gauge invariant functionals, we make use of the concept of gauge 
invariant projection, defined as follows: Let / = L] be an arbitrary functional 
of the phase-space fields. Then its gauge-invariant projection V(I)[ir,L] is defined 
by: 

V(I)[tt,L} = 1 JvuI[^,L], (30) 

where n u is the gauge-transformed of 7r by the gauge group element u>(x) (for exam- 
ple, in 2 + 1 dimensions, 7r^ = %j + ej k D k uj) and the functional integration is over 
all the possible configurations for u. The normalization factor M is just the volume 
of the gauge group: M = J Vu. It is then easy to see that the gauge invariant 
projection of an arbitrary functional is indeed gauge invariant: 

V(I)K,L] = V(I)[tt,L} , (31) 

and that V is a linear projection operator: 

V(X 1 I 1 + X 2 I 2 ) = MVih) + X 2 V{I 2 ) , 

V 2 = V , VJ. (32) 

A functional F is gauge invariant if and only if V(F) = F. This can be shown to 
be equivalent to saying that F belongs to the image of V. We then construct the 
most general gauge-invariant functional by applying V to an arbitrary functional. 

In 2 + 1 dimensions, we further decompose the momentum as 

7Tj(x) = DjCt(x) + e jk D k (3(x) , (33) 

(where a and (3 are scalar and pseudoscalar, respectively) to show that 

V(I)[n,L] = V(I)[a,/3,L] = ±fJvuI[D ja + e jk D k (P + u),L] 

= I[a,0,L] (34) 
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where the last line was obtained by performing the shift lu — > u — j3. ([34]) shows that 
any gauge invariant functional is independent of f3; the reciprocal is immediate. The 
conclusion can be put as follows: The general gauge invariant functional depends 
arbitrarily on L, an on tt only through the combination DjiTj. 

This result is generalizable to 3 + 1 dimensions. F is shown to depend only on 
DjiTj and Lj, by using the same argument as in the 2 + 1 case. The decomposition 
of TTj is now 

TTj(x) = Dj(x(x) + € jk iD k Pi(x) , (35) 

and the j3 dependence is removed as before by a shift in u. The only difference 
appears in the actual construction of the projection operator, which appears to be 
ill-defined at first sight. This is so because the gauge transformations in d = 3: 

71 ~j = n j + € jkl D k 0Ji , (36) 

are invariant under ujj(x) —>■ oJj{x) + Dj\(x), for any A. This produces an infinite 
factor when one integrates over u in the definition ([3(]) of V(I). Of course, this factor 
is also present in A/", but to explicitly cancel them on needs to 'fix the gauge' for 
the integration over u. A convenient way to do that is by using the Faddeev- Popov 
trick, which gives the 'gauge fixed' projector 
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JVudesbM f [u]S[f(u)] 
x J Vudet M f [uW(u}))I[D ja + e ifcl D fc w,] , (37) 



where M/fcu] = j^f(uj x ). We have seen that the gauge invariant functionals depend 
on DjiTj and Lj (the result is indeed true in any number of dimensions). However, 
there is still a degree of redundancy in this description because one is interested 
only in gauge invariant functions on-shell, i.e., on the surface Fjk(L) = 0. Thus we 
do not need the full Lj, but only its restriction to the constraint surface. As it was 
shown in ref. |§, it is possible to solve that kind of equation using a perturbative 
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approach. The main result we need to recall is that that perturbative expansion 
allows one to express L 3 as a function of the scalar d 3 L 3 only. Then we obtain a 
more symmetrical description in terms of the gauge-invariant, scalar variables: 

(D 3 n 3 ) a (x) , (D j L j ) a (x) = djL«{x) . (38) 

Their equations of motion link each other: 

d 

d 

■^(djLj) = -g 2 d j D j (D k 7r k ) 

F jk {L) « 0, (39) 

(where we have included the constraints). They imply the second order equations 

1$ - djDfiDwitk = 

i^t ~ 9jDj)d k L k = 0, (40) 

which show the scalar particle nature of the (only) physical degree of freedom. Let 
us consider in more detail the issue of static solutions in 3 + 1 dimensions. In this 
situation, (39) reduces to 

djLj = , F jk w 
d 3 D 3 {D k K k ) = 0. (41) 

The first two equations in (41) are equivalent to 

L 3 = Ud 3 U ] , djLj = . (42) 

They are exactly the set of equations one gets when considering the Gribov prob- 
lem ]T3| (for the Yang-Mills theory) in the Coulomb gauge, on the orbit of the trivial 
configuration (L 3 =0). It is well known that there are more solutions than just the 
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trivial one, in particular, one obtains the 'fermionic' configurations of the Skyrme 
modef], which verify 

U = I 6jkl tr ^ LkLl ^ = ±2 • ( 43 ) 

Once a particular solution of ([43]) is obtained, it can be inserted in the last equation 
of fl4"l"| ) to get an equation for 7T. Note that the momenta should then satisfy 

Df*i = fo , (44) 

where fo is a zero mode of the operator djDj (of course, the trivial solution %j = 
is included). For each Gribov solution Lj, there will be a non-trivial zero mode for 
this operator, and then a non-zero solution for the momenta. These solutions can 
be compared to the static solutions of the usual non-polynomial formulation. To do 
that we must regain the field Lq, which was eliminated by using the second-class 
constraints. That is very simple, since in fact Lq is equal to a constant times Dj7Tj, 
and then ([44]) implies 

L = fo ■ (45) 

So the family of static solutions in the polynomial version seems to be larger than 
in the usual treatment. Indeed, as Lq = UOqW , a non-zero Lq implies that there is 
a time-dependence for U. Note, however, that such configurations contribute to the 
energy in an amount: 

E(L ) = l -g 2 J 'd'xifo} 2 , (46) 

which is proportional to the norm of the zero mode, and then the minimum en- 
ergy will correspond to the trivial configuration Lq = 0. A simple example of a 
configuration with Lq ^ is: 

Lj(x.,t) = exp(iht)Lj{x) exp(-iht) 

L (t) = h, (47) 

7 The stabilizing term can be added without changing the canonical structure of the model. 
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where h is a hermitian (constant) traceless matrix, and %(x) satisfies (|42]). Thus 
for (|47|), E(L ) = ^g 2 tr(h 2 ) J d 3 x, which is divergent for infinite volume. 

4 Dirac's brackets method 

As an alternative to the previous approach, we apply here the 'Dirac's brackets 
method' to the treatment of the first-class constraints in the 2 + 1 model (it can 
however be straightforwardly generalized to the d + 1 model). It consists in con- 
structing the Dirac's brackets for the set of second-class constraints containing all 
the original first-class constraints plus a suitable set of gauge fixing conditions. We 
choose the canonical gauge fixing functions: 

X a (x) = 7r a 2 (x) = . (48) 

The basic ingredient to calculate the Dirac's brackets is the Poisson bracket between 
X a and G a (x): {x a (x) , G b (y)} = (_Di) ab £(x — y). From this it follows that the only 
non-trivial Dirac's brackets between canonical variables are 

{L a 2 {x)^\{y)} D = (x,a\D^D 2 \y,b) . (49) 

The second one is a complicated non-local function. It is more convenient to take 
advantage of the results of the previous section to work with Lj and Diiri. Then 
the Dirac's brackets become local 

{L?(x),(IW(2/)}d = -D?*(x-y) 

{L^^D^fiynn = -D*6(x - y) . (50) 

5 Conclusions 

The polynomial formulation (0) has an interesting canonical structure. Some of its 
properties are: 
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The system has second-class constraints which can be solved explicitly for some 
coordinates in terms of the others. This leaves the canonical pairs associated 
to the spatial components of a non-Abelian vector field only. 

For d > 1 there remain first-class constraints which form an Abelian algebra. 
They, and the first-class Hamiltonian have essentially the same structure in any 
number of dimensions. However, for d > 2, the constraints are reducible. The 
number of independent constraints is just enough to leave only one physical 
degree of freedom. 

These first-class systems can be regarded as 'duals' of the Yang-Mills model in 
the temporal gauge, in the sense that the constraints in one of the systems 
are non-trivial gauge invariant functions in the other. This duality can be 
generalized to a greater class of first-class systems. 

We also constructed the most general gauge invariant functional explicitly. Note 
that the Gauss-law constraints of the dual Yang-Mills system appear here as 
(non-trivial) gauge invariant objects, verifying the general property discussed 
in the Appendix. 
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Appendix: A duality transformation for first-class 
systems 



The kind of 'duality' that exists between the massless Yang-Mills theory and the 
polynomial version of the non-linear a-model is a particular case of a more general 
concept, which we define in this Appendix. Let us consider a constrained dynamical 
system defined on a phase-space of coordinates qj,Pj, j = 1- ■ ■ N, with first-class 
Hamiltonian H and a complete set of irreducible first-class constraints G a ~ 0, a = 
1, • • • , N. We assume that the first-class constraints satisfy the closed algebra 

{G a ,G b } = g a bc(q,p)G c (51) 

and regarding the Hamiltonian, we impose on it the requirement of having the 
structure 

H = \ F a (q,p)F a (q,p) (52) 
where the functions F a (q,p), a = 1, ■ ■ ■ , N verify the relations 

{F a ,F b } = f abc (q,p)F c 

{G ai F b } = X abc {q,p) F c (53) 

and A is completely antisymmetric with respect to the last two indices. This implies 
that the Poisson bracket of H and each of the G a 's will be strongly equal to zero, 
what is stronger than what we need in a general first-class system. Indeed, Equations 



(52) and (|53|) select among all the possible first-class systems the class which admit 
a duality of the kind we are going to define. 

The associated dual first-class system is defined on the same phase-space, and 
its Hamiltonian and constraints (denoted with a tilde) are defined by 

H = — G a G a 

G a = F a ^0, (54) 
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where F a and G a are the ones introduced in (|5lD , (|5^) and fl5"3|). We then verify 
that the new system is also first-class, since 




g a bc(q,p) G c 



{G a , H} 



V ab (q,p) G b 



(55) 



where: 



9abc{q,p) 



fabc(q,P) 



Vab(q,p) 



Kbc(q,p) G c (q,p) . 



(56) 



Thus evidently this mapping leaves the first-class nature of the system invariant. 
Note however, that the irreducibility of the new constraints is by no means guaran- 
teed. That will depend upon the particular form of the F a 's. An interesting property 
of the new system is that, because of (|53|). 



which proves that the G a 's constitute a set of M independent gauge invariant func- 
tions, which is a very helpful property when one wants to study the classical or 
quantal dynamics of the system. 

Note that the transformation we defined is not necessarily involutive; to guaran- 



The Hamiltonian fl52|) resembles the one of the Yang-Mills system, except for 
the absence of the term quadratic in the canonical momenta. We did not include 
this, neither the corresponding one in the dual, to keep the discussion as general as 
possible. As they are gauge invariant by themselves, their presence or not do not 
alter the essence of the discussion. 

This duality transformation can be interpreted as transforming a gauge-invariant 
theory into another. The unphysical gauge variables of the first theory come to be 



{G a ,G b } « Va,5, 



(57) 



tee that we would need a completely antisymmetric \ a bc in Equation (^3|). 
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the true physical degrees of freedom of the second one, which are actually defined 
on the fibers generated by the gauge group of the first model (note that the second 
theory is not invariant under the gauge group of the first one). Hence it follows 
naturally that the first-class constraints which generate the gauge transformations 
of the first theory are the dynamical variables which describe motions ('translations') 
along the gauge group fiber. 

This duality transformation may be useful in cases where the number of physical 
degrees of freedom is the same in both the original and dual models, as for the Yang- 
Mills and non-linear sigma- model in 2+1 dimensions. Here the physical excitations 
are massless scalar fields for both models. Although it is very difficult to write down 
an explicit dynamics for the physical degree of freedom of the Yang-Mills theory, 
this mapping could make it easier, since the identification of physical variables is 
much simpler in the non-linear sigma-model. 
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